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Abstract 

The equations of motion of two point masses have recently been derived at 
the 3PN approximation of general relativity. From that work we determine 
the location of the innermost circular orbit or ICO, defined by the minimum 
of the binary’s 3PN energy as a function of the orbital frequency for circu¬ 
lar orbits. We find that the post-Newtonian series converges well for equal 
masses. Spin effects appropriate to corotational black-hole binaries are in¬ 
cluded. We compare the result with a recent numerical calculation of the 
ICO in the case of two black holes moving on exactly circular orbits (heli¬ 
cal symmetry). The agreement is remarkably good, indicating that the 3PN 
approximation is adequate to accurately locate the ICO of two black holes 
with comparable masses. This conclusion is reached with the post-Newtonian 
expansion expressed in the standard Taylor form, without using resummation 
techniques such as Pade approximants and/or effective-one-body methods. 
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The aim of this article is to compute the innermost circular orbit (ICO) of point-particle 
binaries in post-Newtonian approximations, and to compare the result with numerical sim¬ 
ulations. For the present purpose, the ICO is defined as the minimum, when it exists, of 
the binary’s energy function E(Q) of circular orbits, where O denotes the orbital frequency. 
The definition is motivated by our comparison with the numerical work, because this is pre¬ 
cisely that minimum which is computed numerically]]. The energy function is given by the 
invariant quantity associated with the conservative part of the post-Newtonian dynamics, 
i.e. ignoring the radiation reaction effects. It can be argued, because the radiation reaction 
damping is neglected, that the importance of the ICO does not lie so much on its strong 
physical significance, but on the fact that it represents a very useful reference point on the 
definition of which the post-Newtonian and numerical methods agree. 

The question of the conservative dynamics of compact binary systems has been resolved 
in recent years at the third post-Newtonian (3PN) approximation, corresponding to the or¬ 
der 1/c 6 beyond the Newtonian force. After the previous work of Refs. [|1].|2[, Jaranowski 
and Schafer [[|, and Damour, Jaranowski and Schafer P|, have applied at 3PN order the 
ADM-Hamiltonian formalism of general relativity. On the other hand, extending the method 
of Ref. [!5[], Blanchet and Faye [|6|,j?f performed a 3PN iteration of the equations of motion 
(instead of a Hamiltonian) in harmonic coordinates. In the latter approaches, the compact 
objects are modelled by point particles, described solely by two mass parameters rri\ and 
m 2 . The end results are physically equivalent in the sense that there exists a unique trans¬ 
formation of the particle’s dynamical variables that changes the 3PN harmonic-coordinates 
Lagrangian of de Andrade, Blanchet and Faye || into another Lagrangian, whose Legendre 
transform is identical with the 3PN ADM-coordinates Hamiltonian of Damour, Jaranowski 
and Schafer J3|. 

Post-Newtonian computations of the motion of point particles face the problem of the reg- 


1 In particular, we do not define the ICO as a point of dynamical (general-relativistic) unstability. 
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ularization of the infinite self-field of the particles. The regularization scheme of Hadamard, 
in “standard” form, was originally adopted in the ADM-Hamiltonian approach |3|]. Then an 
“improved” version of this regularization was defined in Refs. [7J and applied to the com¬ 
putation of the harmonic-coordinates equations of motion |6]. Unfortunately, it has been 
shown that the Hadamard regularization, either in standard or improved form, leaves un¬ 
specified one and only one numerical coefficient in the 3PN equations of motion, o; s t a tic in 
the ADM-Hamiltonian approach, A in the harmonic-coordinates formalism. The parameter 
^static can be seen as due to some “ambiguity” of the standard Hadamard regularization, 
while A appears rather like a parameter of “incompleteness” in the improved version of 
this regularization. However, these constants turned out to be equivalent, in the sense that 

B01 


3 1987 

A UJ of cifip ’ 

11 3080 


( 1 ) 


ft has been argued in Ref. |J that the numerical value of u; s tatic could be ~ —9, because for 
such a value some different “resummation” techniques, when they are implemented at the 
3PN order, give approximately the same numerical result for the ICO. Even more, it was 
suggested |§ that u; sta tic might be precisely equal to u;* tatic = — y + |^7r 2 ~ —9.34 (corre¬ 
sponding to A* ~ 1.90). But, more recently, a computation of a; static has been performed by 
means of a dimensional regularization, instead of the Hadamard regularization, within the 
ADM-Hamiltonian formalism |T0|, with result 


^static = 0 A = -7^ - -0.64 . (2) 

We adopt in this article the latter value as our preferred one, but in fact it is convenient to 
keep the ambiguity parameter unspecified, and to investigate the behaviour of the solutions 
for different values of A or a; sta ti C . For instance, we shall keep an eye on the values u;* tatic — 
—9.34 and also A = 0 u; sta ti C — —2.37. The latter case corresponds to the special 
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instance where certain logarithmic constants associated with the Hadamard regularization 


in harmonic coordinates do not depend on the masses [Q. Notice that the result (^j) is 
quite different from u;* tatic ~ —9.34 : this suggests, according to Ref. ||, that different 


resummation techniques, viz Pade approximants [11] and effective-one-body methods [12 


which are designed to “accelerate” the convergence of the post-Newtonian series, do not in 
fact converge toward the same “exact” solution (or, at least, not as fast as expected). 

Let us now compute the ICO of two point particles (modelling black holes) at the 
3PN order thanks to the previous body of works []l]— . The circular-orbit binding energy 
E (in the center-of-mass frame), and angular-momentum J, are deduced either from the 
3PN harmonic-coordinates Lagrangian || or, equivalently, from the 3PN ADM-coordinates 
Hamiltonian [f|| (we neglect the 2.5PN radiation damping). These functions are expressed in 
invariant form (the same in different coordinate systems), i.e. with the help of the angular 
orbital frequency H. The 3PN energy (per unit of total mass M), describing “irrotational” 
circular-orbit binaries, is 
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All over this paper we pose G — c — 1. Mass parameters are M = mi + m 2 , and the 
symmetric mass ratio v = mpm^/M 2 such that 0 < v < 4 , with v — | in the equal-mass 
case and v —* 0 in the test-mass limit for one of the bodies. The 3PN angular momentum, 
scaled by M 2 , reads 
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The variations of the energy and angular momentum of the binary on the circular orbit 
during the inspiral phase obey the evolutionary (or “thermodynamic”) law 
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which is equivalent, via the energy and angular-momentum balance equations, to the same 
relation but between the corresponding gravitational-wave fluxes at infinity. From Eq. (i. 
we see that the points of extremum for E and J are the same. In the limit v —> 0, Eqs. 
(l)-i) reduce to the 3PN approximations of the known energy and angular momentum of 
a test particle in the Schwarzschild background : 

ALL = „{(!- 2(M(l) 2 / 3 ) (1 - 3(Mfl ) 2 ' 3 ) _1/2 - 1} , (6a) 

ALL = (1 - 3(MSi) 2/s y 1/2 . (6b) 

We recall that in this case the location of the ICO is given by MOj^q = 6 -3//2 , with Ef$Q = 
uM — lj and = ^M 2 a/T2. 

The straightforward post-Newtonian method we follow in this article can be justified by 
the following arguments. At the location of the ICO we shall find that MOjco is of the 
order of 10%. Therefore, we expect that the 1PN approximation will grossly correspond to 
a relative modification of the binding energy of the order of v 2 ~ (MOico) 2 ^ 3 i-e. 20%; and 
similarly that the 2PN and 3PN approximations will yield some effects of magnitude about 
5% and 1% respectively. Consequently the post-Newtonian method should be adequate in 
the regime of the ICO, provided that it is implemented up to the 3PN order, so as to be 
accurate enough. On the other hand, we see that the 1PN order should yield a rather poor 
estimate of the position of the ICO. 

Let us now confirm these estimates with the numerical values for the post-Newtonian 
coefficients in the energy function (j3j) . As we see from TABLE [I|, in the case of comparable 
masses and of our preferred value (^j) for the ambiguity parameter, the absolute values of 
the post-Newtouian coefficients are roughly of the order of one (they do not apparently 
increase with the order of approximation). This means that the previous estimates are 
essentially correct. In particular the 3PN approximation should be close to the “exact” 
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value for the ICO. The post-Newtonian series seems to “converge well” (in the case where 
v — | and o; stat i C = 0), with a “convergence radius” of the order of one, i.e. at a much higher 
frequency than the frequency of the ICOjj. By contrast, we recover in TABLE I the well- 
known result (see e.g. |1^,[L4|]) that in the perturbative case v —> 0 the post-Newtonian series 
converges slowly : the coefficients increase roughly by a factor 3 at each post-Newtonian 
order, reflecting the fact that the radius of convergence of the series is This is clear from 


the exact expression (|6a|) , in which the pole at the value | corresponds to the light ring 
of the Schwarzschild metric. Thus the post-Newtonian method is not very appropriate to 
the case v = 0, where even the 3PN order would rather poorly approximate the ICO. The 
situation is therefore the following : in the case of comparable masses, we do not have the 
exact solution, but fortunately the straightforward post-Newtonian approach is expected to 
be accurate; in the perturbative limit v = 0, the post-Newtonian series is poorly convergent, 
but gladly this does not matter because we know the exact results (||). 

Having thus justified the validity of our approximation, we look for the point at which 
both E(Q) and J(Q) take some minimal values E ICO = E(ff ICO ) and Jico = ^(^ico)- As we 
see from Eq. (j3j), at the 3PN order E(Q) is a polynomial of the fourth degree in the frequency- 
parameter x = (MO) 1 2 ^ 3 . Therefore, the value of the minimum, xico = (MJlico) 2 ^ 3 , must 
be a real positive solution of an algebraic equation of the third degree (in general) : 


1 + ax + (3x 2 + 'yx 3 = 0 . (7) 

The coefficients are straightforwardly obtained from Eq. (|3|) as 

2 Actually the post-Newtonian series could be only asymptotic (hence divergent), but neverthe¬ 
less it should give good results provided that the series is truncated near some optimal order of 

approximation. In this article we assume that 3PN is not too far from that optimum. 
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The regularization constant A enters only the third-degree monomial (3PN order). Let us 
describe, in a qualitative way, the existence of solutions of Eq. (^). We find that the 
equation does not always admit a unique real positive solution, nor even several of them. 
This depends, for a given choice of the mass ratio z/, on the constant A. When A happens to be 
smaller that some “critical” value Ao(zz), depending on z/, there is no (real positive) solution, 
and therefore there is no ICO at the 3PN order. When A is between Ao(zz) and another 
“critical” value Ai(z/), also depending on z/, we obtain two real positive solutions. In this 
case, the energy function admits two extrema, a minimum and a maximum. The maximum 
occurs at a higher frequency than the minimum of the ICO, and is to be discarded on 
physical grounds (the corresponding frequency is generally too high, e.g. higher than M~ 1 , 
for being of physical interest). Finally, when A is larger than Ai(z^), there is one and only 
one real positive solution : aqcch and this is a minimum of the energy. The latter regime, 
where the circular-orbit energy admits a unique extremum, which is a minimum (like for the 
Schwarzschild metric), is the simplest on the physical point of view. The interesting values 
of A are located in the regime where A > Ai(z/) (for irrotational binaries). We summarize 
our discussion in FIG. [Tj. 

It is not difficult to determine analytically the functions A 0 (zz) and Ai(z/). Indeed, A 0 (^) 
represents simply the minimal value of the function xico — 1> A(zy xico) (see FIG. |l!). Using 
also Eq. ( 0 ), we readily find the mathematical relation defining Aq(zz) : 


A — A 0 (z/) 


7(UA) = - 


[a 2 (u) - 3/3(z /)] 3/2 - a 3 {v) + ^a(z/)/3(z/) 


( 9 ) 


from which the explicit expression of Ao(zz) can be found using Eqs. (] 8 j). On the other hand, 
the function Ai(z/) is determined by the cancellation of the third-degree coefficient in the 
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equation (|7|), i.e. 


A = Ai(z/) 


7 (z/, A) = 0 . 


( 10 ) 


The expression of Ai(z/) then follows from using Eq. (|5c|). For allowed values of v G 
]0, |], we find that both Ao(u) and Ai(z/) are increasing functions of u, with maximal values 
Ao(^) — —2.2 and Ai(^) ~ —0.96, and satisfy Ao(u) —»► —oo and \\[y) —> —cx) when v —> 0. 
Furthermore we always have Ao(u) < Ai(u). This analysis shows that in the case of our 
preferred value A = «|£|Z ~ —0.64, as well as in the cases where a; sta tic = —9.34 and A = 0, 
the energy function E(Q) given by Eq. ([|), for any mass ratio u, admits a unique extremum, 
which is a minimum, at some flico (f° r corotating binaries we shall find a minimum and 
also a maximum at very high frequency). We show in FIG. the graph of E(£2) for equal 
masses and a; sta tic = 0. Anticipating on our discussion below, it is interesting to compare 
FIG. |2] with the result of the numerical simulation provided by the figure 16 in Ref. El- 

In TABLE |I]we present the values of the calculated frequency Ojccu the corresponding 
energy E\co and angular momentum JicO; at the 1PN and 2PN orders, and at the 3PN 
order in the three cases where u; static = 0, A = 0, and u; sta ti C = —9.34. The 1PN and 2PN 
approximations are defined by the obvious truncation of Eqs. (|3])-([4|)- Notice how close 
together already are the 2PN and 3PN approximations (however, the 1PN order seems to be 
quite inadequate). Let us now show that the 3PN approximation, in standard form (Taylor 
approximants), appears to be very good to locate the turning point of the ICO, in the sense 
that the prediction for that point is close to the recent result of numerical relativity. 

A novel approach to the problem of the numerical computation of binary black 
holes in the pre-coalescence stage, has been proposed and implemented by Gourgoulhon, 


Grandclement and Bonazzola [pLo| , [TG|l . This approach uses multi-domain spectral methods 
T7|] , and is based on two approximations, the first one is essentially “technical”, the other 


one is “physical”. The technical assumption (which could be relaxed in future work) is the 
conformal flatness of the spatial metric : 7 ^ = T 4 ^-. On the other hand, an imposed “he- 
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lical” symmetry constitutes an important physical restriction to binary systems moving on 
exactly circular orbits. By helical symmetry we mean that the space-time is endowed with 
a Killing vector held of the type ^ where and ^ denote respectively the 

time-like and space-like vectors that coincide asymptotically with the coordinate vectors of 
an asymptotically inertial observer. A crucial advantage of the helical symmetry, especially 
in view of the comparison we want to make with the post-Newtonian calculation, is that 
the orbital frequency 0 is unambiguously defined as the rotation rate of the Killing vector. 
Thanks to these approximations, Gourgoulhon, Grandclcment and Bonazzola []15|,|I(]] were 
able to obtain numerically the energy and angular-momentum along the binary’s evolution¬ 
ary sequence, i.e. maintaining Eq. (|]) along the sequence, and to determine the minimum 
of these functions or ICO. 

The numerical calculation reported in Refs. (l^JT^| has been performed in the case of coro- 
tating black holes, which are spinning with the orbital angular velocity 0. We must therefore 
include within our post-Newtonian treatment the effect of spins]], appropriate to two Kerr 
black holes rotating at the orbital rate 0. By combining the formula of Christodoulou and 
Rufhni : m 2 = mf rr + -/4-, with the known relation between the black-hole spin and its 


angular velocity]] : S = 2m 3 4 Q 


1 + 


s 2 


, we obtain the total mass m and spin S of 


each of the corotating black holes in terms of their irreducible mass m nr . 


m h 


m = 


S = 


a/1 - 4(m irr O) 2 
4mf rr f2 

a/ 1 - 4(m irr 0) 2 


~ m irr + 2 m/O 2 , 


~ 4m/.Q 


( 11 a) 

(lib) 


3 The importance of the effect of spins in corotating systems of neutron stars, for which the ICO is 
usually determined by the hydro dynamical instability rather than by the effect of general relativity, 
is well known 

4 More precisely the angular velocity is defined as the one of the outgoing photons that remain for 


ever at the location of the horizon; see Eq. (33.42b) in Ref. [19]. 
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The irreducible masses are precisely the ones which are held constant along the evolutionary 


sequences calculated numerically in Refs. [|15|,|Kj. Therefore our first task is to replace all 
the masses parametrizing the sum M + E, where M = mi + m .2 is the total rest mass-energy 
and E is the 3PN binding energy given by Eq. (]3]), by their equivalent expressions, following 
Eq. (|lla| ), in terms of the two irreducible masses. It is clear that the leading contribution 
is that of the kinetic energy of the spins and will come from the replacement of the rest 
mass-energy M; from Eq. (|lla|) we see that this effect will be of order ff 2 in the case of 
corotating binaries, which means by comparison with Eq. (||) that it is equivalent to an 
“orbital” effect at the 2PN order, ffigher-order corrections in Eq. ( |lla|) will behave at 
least like f2 4 and correspond to the 5PN order at least, negligible for the present purpose. 
In addition there will be a subdominant contribution, of order ff 8 / 3 or 3PN, coming from 
the replacement of the masses into the “Newtonian” part, oc f2 2//3 , of the binding energy E 
[see Eq. (|3])]. At the 3PN approximation we do not need to replace the masses into the 
post-Newtonian corrections in E. Our second task is to include the relativistic spin-orbit 
(S.O.) interaction. In the case of spins Si and A 2 aligned parallel to the orbital angular 
momentum (and right-handed with respect to the sense of motion) the S.O. contribution to 
the energy reads 00 


E s _ o. = -vM (MO ) 5/3 


4 rn\ 
3 M 2 


Si 


+ V - o + X 


m; 
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3 M 2 


n i S * 

+ u I -2 


( 12 ) 


As can immediately be infered from S ~ 4m 3 0, which is deduced from Eq. (|llb| )p|, in the case 
of corotating black-holes the S.O. effect is of order 3PN and therefore must be retained at the 


present accuracy [with this approximation, the masses in Eq. (|T2D can be chosen to be the 
irreducible ones]. By contrast, the spin-spin (S.S.) interaction turns out to be much smaller, 
equivalent to the 5PN order for corotating systems. Considering all the contributions present 


5 The moment of inertia of the Kerr black hole in the limit of slow rotations is I = 4m 3 , in 


accordance with Eq. (2.61) in Ref. [22 
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with the 3PN accuracy, we thus obtain three terms : (2 — 6u) (Mfl) 2 coming from the kinetic 
energy of the corotating spins; (—|z/ + iA) (MQ) 8 ^ 3 due to a coupling between the spin 
kinetic energy and the orbital energy; (— yz/ + 12z/ 2 ) (Mf2 ) 8 ^ 3 due to the S.O. interaction 
(0)- Numerically the kinetic energy of the spins will dominate the other effects. Hence the 
supplementary energy that is due specifically to the corotation reads 


E COTOt (tt) 

M 


= (2 - 6 v) (MO ) 2 + + 13z/2 ) ( Mfi ) 8/3 • ( 13 ) 


The total binding energy of the corotating binary is the sum of Eqs. (|3[) and ( 0 ). Notice 
that we must now understand all the masses in Eqs. (|3]) and (jl3[) as being the irreducible 
masses (we no longer indicate the superscripts “irr”), which stay constant when the binary 
evolves following Eq. (|5|). 

In TABLE [Hi] we present our results for Eico and Qico of a corotational binary. Since 
E corot , given by Eq. ([13]), is at least of order 2PN, the result for lPN corot is the same as for 


1PN in the irrotational case; then, obviously, 2PN corot takes into account only the leading 
2PN corotation effect (i.e. the kinetic energy of the spins), while 3PN corot involves also, in 
particular, the corotational S.O. coupling at 3PN order. In FIG. ^jwe plot Eico versus Oico, 
computed with and without the corotation effect, and compare the values with the result 


obtained by numerical relativity under the assumption of helical symmetry |16|. As we can 
see the 3PN points, and even the 2PN ones, are rather close to the numerical value. As 
expected, the best agreement is for the 3PN approximation and in the case of corotation'] : 
i.e. the point 3PN corot . However, the 1PN approximation is clearly not precise enough, but 
this is not very surprising in this highly relativistic regime where the orbital velocity reaches 
v ~ (MLico ) 1 / 3 ~ 0.5. Summarizing, we find that the location of the ICO computed by 


3 We have checked that our best value, given by 3PN cor , is not significantly modified numerically 


when we add the higher-order spin effects in Eq. (13) up to the 5PN order, i.e. including, in 
particular, the S.S. interaction. 
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numerical relativity, under the helical-symmetry approximation, is in good agreement with 


post-Newtonian predictions. This was already pointed out in Ref. (ldj from the comparison 
with Pade and EOB methods. This constitutes an appreciable improvement of the previous 
situation, because we recall that the earlier estimates of the ICO in post-Newtonian theory : 
MO IC o — 0.06 and Ei CO /M ~ —0.009 PB[ , and numerical relativity : Mfl ICO — 0.17 and 
Eico/M ~ —0.024 


5[, strongly disagree with each other, and do not match with the 
present 3PN results (see Ref. [[R| for further discussion). 

Let us emphasize that our computation has been based on the standard post-Newtonian 
approximation, expanded in the usual way as a Taylor series in the frequency-related param¬ 
eter x = (Mf2) 2 ^ 3 [see Eqs. ©-(§) and (0)], without using any resummation techniques. 
In FIGS. f|-|]we display our Taylor-series-based values for Eico and Jico (they are indicated 
by the marks 2PN and 3PN), and contrast them with some results obtained by means of 
resummation techniques at the 3PN order : Pade approximants [|TT|,|9| and effective-one-body 
(EOB) methods [fT2|.|0|| . All these results agree rather well with each other, and, as we have 
seen, even the 2PN (Taylor) approximation does well. 

A point we make is that the sophisticated Pade approximants give about the same results 
as the standard post-Newtonian expansion, based on the much simpler Taylor approximants : 
indeed, see in FIGS. Q-[| the points referred to as the e and )-met hods, which are 3PN Pade 
resummations built respectively on the energy and angular-momentum |J. For the case 
at hands - equal-mass binaries —, there is apparently no improvement from using Pade 
approximants. Nevertheless, it is true that in the test-mass limit u —> 0 the Pade series 
converges rapidly toward the exact result |TT[. For instance, the Pade constructed in this 
case from the 2PN approximation of the energy already coincides with the exact expression 
for the Schwarzschild metric [given by Eq. (|6a|)]. But, the results of FIGS. Q-|5] suggest that 
this interesting feature of the Pade approximants is lost when we turn on v and consider the 
equal-mass case v — Notice also that the 2PN versions of these Pade, which are given in 
the table I of Ref. differ much more significantly from the corresponding 3PN ones than 
in the case of Taylor. For instance, the 2PN e-method yields the values Mflico — 0.09 and 
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Eico/M ~ —0.016, which respectively differ by about 36% and 22% with the frequency and 
energy given by the e-method at 3PN. In the case of Taylor, the same figures are only 6% 
and 3%. Thus, on the point of view of the “Cauchy criterium”[]. the Taylor series seems to 
converge better that the Pade approximants (for equal masses). 

It is a pleasure to thank Eric Gourgoulhon for informative discussions, Alessandra Buo- 
nanno and Gillcs Esposito-Farese for useful remarks, and a referee for valuable comments. 


'The Cauchy criterium for the series a n is the fact that \a n — a m \ —> 0 for any n and m. 
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TABLES 



Newtonian 

1PN 

2PN 

3PN 

^ — 4 ^static — 0 

1 

-0.77 

-2.78 

-0.97 

V = 0 

1 

-0.75 

-3.37 

-10.55 


TABLE I. Numerical values of the sequence of coefficients of the post-Newtonian series com¬ 
posing the energy function 




Afflico 

E\co 

M 

Ago 

~w 

1PN 


0.522 

-0.0405 

0.621 

2PN 


0.137 

-0.0199 

0.779 

3PN 

^static — 0 

0.129 

-0.0193 

0.786 

3PN 

A = 0 

0.116 

-0.0184 

0.798 

3PN 

^static = 9.34 

0.095 

-0.0166 

0.824 


TABLE II. 

Parameters for the ICO of equal-mass {v 

= \) binary systems. 




Mflico 

Aioo 

M 

lPN corot 

0.522 

-0.0405 

2PN corot 

0.081 

-0.0145 

3PN corot o; static = 0 

0.091 

-0.0153 


TABLE III. Parameters for the ICO of corotational equal-mass binary systems. 
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FIGURES 


X ico 



FIG. 1. The possible solutions as a function of the regularization constant A. There is no 
solution when A < Ao(z'), two possible solutions when Ao(z') < A < Ai(zz) [which become degenerate 
at A = Ao(^)], and a unique solution when Ai(z/) < A. The upper branch, existing between Ao(t') 
and the vertical asymptote at A = Ai(z/), is actually a maximum of the energy. 
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FIG. 3. Results for Efco versus Rico i n the equal-mass case. The asterisk marks the result 
calculated by numerical relativity. The points indicated by 1PN, 2PN and 3PN are computed from 
Eq. ©, and correspond to irrotational binaries. The points denoted by lPN corot , 2PN corot and 
3PN corot come from the sum of Eqs. (&) and (13), and describe corotational binaries. Both 3PN 
are 3PN corot are shown for astatic = 0. 
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FIG. 4. Results for Efco i n terms of fhco in the equal-mass case. The e and j-methods are 
Pade approximants at the 3PN order. EOB refers to the effective-one-body approach at the 3PN 
order. The points marked by 2PN and 3PN correspond to the standard Taylor post-Newtonian 
series (this work). The results for Pade, EOB and Taylor are for irrotational binaries. 
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FIG. 5. Same as FIG. [| but for the angular-momentum Jico- 
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